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Abstract 



We study the effects of correlations on a one dimensional ring threaded by a 
uniform magnetic flux. In order to describe the interaction between particles, 
we work in the framework of the U oo Hubbard and t-J models. We focus 
on the dilute limit. Our results suggest the posibility that the persistent 
current has an anomalous periodicity 4>o/p, where p is an integer in the range 
2 < p < N e (N e is the number of particles in the ring and 4>o is the flux 
quantum). We found that this result depends neither on disorder nor on the 
detailed form of the interaction, while remains the on site infinite repulsion. 
PACS number(s): 05.30.Fk, 73.20.Dx, 71.30.+h 

I. INTRODUCTION 

Persistent currents in metal rings enclosing a magnetic flux were first studied in the 
sixties ]]]]. In ideal systems the periodicity of the magnetoconductivity is given by the flux 
quantum O — hc/e. In 1981, Al'tshuler, Aronov, and Spivak |2j renewed the interest in the 
topic predicting that in highly disordered systems the period of the magnetoconductivity 
was 0o/2. This effect was then confirmed experimentally by Sharvin and Sharvin ||. Some 
years later, Biittiker and co-workers predicted persistent currents in a one dimensional loop 
of normal metal driven by an external time-dependent magnetic flux with elastic ||J] and 
inelastic scattering. 

The predicted current in multichannel normal metal rings was analyzed as a function of 
disorder and temperature || and it was found that disorder strongly reduces the amplitude 
of the persistent current. 

In 1990, Levy, Dolan, Dunsmuir, and Bouchiat |7j found experimental evidence of these 
currents. In this experiment, the response of an ensemble of mesoscopic copper rings was 
measured as a function of the enclosed magnetic flux. They observed a periodicity of half a 
flux quantum. This result motivated a great activity concerning the problem of persistent 
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currents in an ensemble of rings, including the interaction between electrons || and disorder 
PJT0|j. In 1991, Chandrasekhar et al. measured the current in a single, isolated gold 



ring. In Refs. J7| and |nj the rings were metallic, but in 1993 Mailly et al. [12j reported 
measurement of the persistent current in a semiconductor ring. Recently, new experiments 
were made in those rings [DJ emphasizing the differences between isolated and connected 



geometries. 

For a metallic loop with impurities and in the difussive regime, the magnitude of the 
current is expected to be ~ (evF/L)rj/L || (where L is the perimeter of the ring, vf is the 
Fermi velocity and rj is the elastic mean free path). This is observed in ensembles of metallic 
rings |7j and in semiconductor rings [|12J but in an isolated ring a persistent current of 



one to two orders of magnitude larger than that predicted by theory was measured. This 



fact has stimulated many recent theoretical works with controversial results [^4| -[T9[ . 

The aim of this work is to present results in the limit of strong coupling and low carrier 
density which are independent of the properties of the material. The paper is organized as 



follows. In Sec. [TJ we present the model and discuss its different limits. In Sec. jTTTj, we show 



the corresponding results and finally, in Sec. [TV], we give their physical interpretation and 
the conclusions. 

II. MODEL 

We study a strictly one dimensional ring threaded by a magnetic field. This is such 
that produces a flux <p concatenated by the ring. In this system the ground state carries 
a steady current, which is periodic in the magnetic flux threading the loop, with period 
0o — he/c, the flux quantum. The current arises from the boundary conditions imposed by 
the magnetic flux. 

To study the properties of a system of interacting electrons in a ring with N s sites we 



use the model proposed by Hubbard [20]. The Hamiltonian of the extended Hubbard model 



is: 
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N s -1 



Hu = -t( E iW+ E e^^c^J+h.c. 



n 1 , ct= ^ 

AT S N s N a r 

+ E e i (^>t + ™a) + ^ E ™*>T ^.i + E E ^™ ^i,afli + m,a' (1) 
i=l i=l i=l m=l 

o-,o-'=T,| 

where Cj iCT (cJ CT ) are the annihilation (creation) operators which annihilate (create) an electron 
of spin a on a site i of the ring, n i)CT = c\ a Ci )C7 is the number operator that counts the number 
of electrons of spin a on the site i, and N s is the number of sites in the chain. 

The first term represents the kinetic energy describing the hopping of an electron from 
one site to a nearest neighbour site with hopping matrix element t. We shall put t — 1 
throughout this paper, thus fixing the energy scale. 

The energies q represent the disorder in the ring, and they can be any number in the 
range between [— W/2, W/2] with equal probability. 

The last term allows us to represent a long range contribution (r gives the extent) where 
V m is the strength of the interation. When this term is zero, we have the Hubbard model in 
which the long range contribution due to the Coulomb interaction is supposed to be screened 
and only is retained when the electrons are in the same site giving an additional energy U 
when the site is doubly occupied. 

Here we will consider neither the interaction through phonons nor any other solid exci- 
tation and therefore we restrict ourselves to the case U > 0. 

If we are in the particular case where {7 = and W = 0, the ground state consists in 
doubly occupied levels filling up to the Fermi level. In this state the total spin is the lowest 
possible. When {7^0 the Coulomb repulsion tends to reduce double occupancy, increasing 
the total spin of the ground state. In this way, the Coulomb energy is reduced in a quantity 
proportional to U but the kinetic energy is enhanced. Then, it is found that the ground 
state is determined by a competition between the Coulomb and kinetic energies. In one 
dimension and with adequate boundary conditions it is possible to demonstrate analitically 
21]] that the ground state has always the minimum total spin. In our case, where the system 



is finite and with periodic boundary conditions, we will see that the result depends on the 



number of particles in the ring. 

On the other hand we consider the extreme case in which the interaction intrasite is 
infinite. The U — oo limit is obtained by projecting onto a subspace where the double 
occupancy is not allowed. The corresponding projection operator P is: 

p = n (i - ( 2 ) 

8=1 

Therefore, the Hamiltonian of the one dimensional U = oo Hubbard model is: 

N s 

Hoc = -tP £ c{ a c l+ i,.P (3) 

a,i=l 

This restriction drastically reduces the size of the Hilbert space explored from a given 
Ansatz by the action of Hoo . In one dimension there is another restriction on the size of the 
available Hilbert space. It comes from the fact that the relative positions of particles having 
different spin must be conserved (up to cyclic permutations) under the action of H^. 

Another possibility to study the large U limit is given by the t-J model which allows us 
to work within the same Hilbert subspace of the U = oo case ||22|| . The Hamiltonian of the 
t-J model is a particular case of the t-J t -J z model 

N s Ns 

H t ^ Jt . Jz = H^+J t £ S+Sr +1 +J Z ^ S?S* +1 +c.c (4) 



( a ) 1 (6) 1 (c) 



when J t = J z = %jr. 



The terms (b) and (c) have a different effect when they act on a state: the (b) term 
connects states having a different spin order on the chain, and the (c) term, being diagonal, 
connects the state with itself. They will be discussed in the next sections. In this way, 
we can study the problem for large U perturbatively beginning from the U = oo Hubbard 
model. 



III. RESULTS 



In this section we show exact results obtained with the Lanczos algorithm. The problem 
of the Lanczos method is that in its implementation it is necessary to deal with vectors of 
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length equal to the size of the Hilbert space under consideration. Since this size usually 
increases exponentially with the number of sites of the lattice, the technique is restricted 
to small chains. This method allows us to obtain the ground state and to compute the 
corresponding current. In this paper we shall concentrate on ground state properties only, 
so we will assume that the energy level separation is in general much larger than the thermal 
energy. In this case the current can be obtained through j = —dE/d(f>. 

A. Hubbard model with finite U 

The persistent current of the non-interacting ring can be understood if we assume that the 
one particle energy levels move along the free particle dispersion relation as the magnetic 
field increases, reducing the energy difference between the levels that correspond to the 
wave vectors —i/N s and (i — 1)/N S . Therefore, if the system possesses 2n particles there 
is an accidental degeneracy between the states with S = and S — 1 (S is the total 
spin of the many-body state). This occurs at a magnetic flux <p c = (fio/2, in the case 
where n is an odd integer or at <p c = when n is an even integer. The current behaves 
monotonically as a function of the flux, having a discontinuity at C . The Fourier spectrum 
of this function has components which correspond to periodicities (fio/m (where m is an 
integer), the most important of them being at m — 1. The Coulomb interaction shifts these 
accidental degeneracies to other values of C : < <p c < (fio/2. At this point there is a 
transition to an S — 1 ground state, because the Coulomb energy gained is greater than 
the kinetic energy lost by increasing the total spin. This transition can also be seen in the 
energy curve as a sharp peak at the transition point. 

If U is large enough, transitions to states with S > 1 are also possible. In any case, the 
number of peaks in the energy curve is increased; then, the Fourier components with m > 1 
are enhanced and the fundamental periodicity is smaller than that of the free particle case. 
This is shown in Fig. I. 

If the number of particles in the ring is an odd number and the flux is close to the 
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accidental degeneracies, there is no change of the total spin. This is so because a change 
of the total spin would imply filling levels of higher kinetic energy than in the even particle 
case (due to the Pauli exclusion principle). So, at least for the smaller values of U, there are 
no peaks in the energy curve. Nevertheless, the transitions will occur at large values of U, 
large enough so as to compensate the kinetic energy increase. This is shown in Fig. 2 where 
we can see that the peaks appear at values of U greater than those of the even particle case 
(see Fig. 1). 

The effects of disorder on these properties are two: it breaks translation symmetry and 
favours states with minimum total spin (the lower energy sites on the chain tend to be 
doubly occupied). So, in this case there are neither accidental degeneracies nor transitions, 
and the fundamental periodicity remains to be O (these conclusions are valid in the weak 
interaction regime). 

In order to study the limit of large U, a natural way is to consider the situation in which 
U = oo, as we will see in the next section. 

B. The U = oo Hubbard model 

Here, we review the U = oo Hubbard model which provides a natural scenario to study 
the t-J model. 

It is well known that infinite one dimensional systems present charge and spin separation. 
In other words we can say that the movement of holes leaves the relative spin order in the 
chain unaltered. Then, the Hamiltonian has a block form, each block corresponding to a 
definite spin order. However, when the system is finite and closed, a particle going from 
the last to the first site of the loop modifies the spin arrangement by a cyclic permutation. 
Even in that case the Hamiltonian has a block form, each block corresponding to a subspace 
where the spin order of the states can only differ by cyclic permutations of the particles. 
These subspaces can be labeled by an integer number F, which is the minimum number of 
cyclic permutations that must be done in order to reobtain the initial state. For example, 
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when we have four particles there are only two subspaces. They are: 



f 1 = 4 — {|nU>;ITUT>;IUTT>;Um>} 



{ltm);imt)} 



(empty sites are irrelevant for this analysis). Clearly, 



<F i \H 00 \F j >=Q if Fi ^ Fj. 



(5) 



For a general number of particles N e , it is possible to show that 2 < Fi < N e . In 
particular if N e is an odd number, the only possibility is F = N e . 

Within each of these subspaces two states that differ by a cyclic permutation are con- 
nected by the boundary terms of the Hamiltonian : 

H b = -te^f^c^ (6) 

In the general case, the Hamiltonian within each F-subspace (H F ) can be written as: 

I fjF,l f e (-i2ir<t>/<p ) . . . f e (-i27r<t>/<Po) \ 

f e (i2n4>/4>o) JjF,2 ... Q 



H 1 



^j> e (i2it<t>/4> ) 







f e (i2n<f>/<p ) jjF,F 



(7) 



Each matrix H F, i is equivalent to a Hamiltonian of N e spinless fermions on a chain with 
open boundary conditions associated with eigenvalues Eq and eingenvectors \j,N e ,a >. 
In terms of this basis of eigenstates the Hamiltonian reads: 

H F = £ E% \j, N e , a) (j, N e , a\+Y, U \h N e , «) (j + 1, N e , (3\ (8) 

3, a j,a,P 

Each term in the first sum describes the matrix H F, K The second sum describes the 
matrix T e ( i27r 4>/M ^ w here t' af3 is given by: 



t' a p = (j,N e ,a\H b \j + l,N e ,(3) 



(9) 



Because of the structure of the if F,J matrices, t' a/3 does not depend on the indices j and 
F. Now, by using 

\j,N e ,a) = -$= y £ e<Pj \P> N e><x) ( 10 ) 
yr p 

where p = (27r/F)n and n is an integer (0 < n < F), we obtain that the Hamiltonian can 
be written as: 

^ = EfEKb^ e ,a)(p,iV e ,«| + cosfp + 27r^)^2t^|p,^ e ,«)(p,iV e ,/5|) (11) 

p 

Three conclusions can be drawn from this formula: 

i) The energy of any eigenstate of H p can be written as: 

E + t'cos(p + 2ir^) (12) 

V 00/ 

with E and t' appropiate constants related to the first and second terms of (11). 

ii) As the flux is increased, there is crossing between the many-body levels that correspond 
to different values of p, as it happens in the free particle case. In addition, the largest number 
of crossing occurs for the largest F value. Therefore, we can say that the ground state always 
belongs to the subspace with maximum F, though possibly degenerate with others of lower 
F . This result is ilustred in Fig. 3 with a numerical example. 

iii) From (|7]) we see that within every F-subspace the problem can be mapped onto a 
tight binding chain with a magnetic flux Ftp. Therefore, the periodicity of the persistent 
current is </>o/F. In particular, the ground state periodicity is (po/N e as it was already 
stated in But we are showing that anomalous flux quantization occurs also in the 
excitated subespaces of the problem, corresponding to the periodicity <po/F (remember that 
2 < F < N e ). Note that in the case where we have an odd number of particles the only 
possibility is F = N e which leads to a persistent current having a periodicity <j)o/N e in each 
subspace. 
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We would like to point out that our analysis is completely general and cannot be modified 
by including diagonal disorder or long range interactions( while the double occupancy of the 
sites remains forbidden), which only affect the coefficients Eq and t a p. In each F-subspace, 
these coefficients are renormalized in the same way (if the interaction does not depend on 
the spin of the particles), so they cannot change the crossings between many body levels, 
driven by the magnetic flux, and the fundamental periodicity remains unaltered. 

C. The t-J model 

Here we consider the effect of a large but finite U by means of the t-J model. 

As a previous step we analyze the t-J z model, obtained by setting J t = in Eq. @. The 
Hamiltonian of this model does not connect states with different F values. We have seen that 
the (a) term favours the maximum value for F, which corresponds to the minimum number of 
antiferromagnetic links. The (c) term favours the minimum value for F, which corresponds 
to the maximum number of antiferromagnetic links. This competition between (c) and 
(a) selects the F-subspace where the ground state belongs, modifying the periodicity of the 
energy and the persistent current, as we show in Fig. 4. Note that the periodicity of a half of 
quantum flux is readily stated as the antiferromagnetic coupling (AFC) is raised, suggesting 
the posibility of superconducting correlations in the model. The effect associated with the 
AFC is to renormalize the eigenvalues Eq and the t aj @ coefficients. This renormalization is 
different within eack F-subspace (because the different number of antiferromagnetic links in 
each one of them), modifying the position of the crossings. The most important armonic of 
the energy curve as function of flux is defined by the lower value of the F number present 
in the curve (see Fig. 4). 

When we include the (b) term, the Hamiltonian connects states that belong to different 
F-subspaces. The effect is to change the spin arrangement. J t << 1 means that the charac- 
teristic time for a change in the spin arrangement is much greater than that corresponding 
to one hole hoppping, so the charge can hop while the spin arrangement remains almost 
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unaltered. Therefore in this limit the periodicity must be the same as in the t-J z model. 
The only effect in this case is the smoothing of the current discontinuities. 

For a strong enough J t , the spin arrangement changes very quickly; then, when the 
charge hops, the phase coherence existing within each F-subspace dissapears. In this case, 
the system recovers a periodicity of one quantum flux. This is shown in Fig. 5. Note that 
even in this case there is a wide region of values of the AFC where the periodicity correspond 
to a half of quantum flux. 



IV. CONCLUSIONS 

In the present work we have studied the effect of strong interactions on the properties of 
mesoscopic rings threaded by a magnetic flux. 

In these systems we predict a fractional Aharanov-Bohm effect with periodicity (fro/p, 
where p is an integer in the range 2 < p < N e . As the strength of the AFC is increased, the 
fundamental periodicity of the persistent current evolves from (fro to (fro/N e . If the particle 
number is even, during this evolution, the system takes the fractional periodicities (fro /p. 
On the other hand, if the particle number is odd, these intermediate transitions do not take 
place, as in this case there is only one value of F, F — N e . To the best of our knowledge, this 
behavior has not been stated before. On a large range of values of the AFC the periodicity 
of the persistent current correspond to a half of quantum flux, suggesting the posibility of 
superconducting correlations in the model. 

Remarkably, our results are independent of the disorder and the detailed form ofthe 
interaction. 

To test our predictions it is necessary to set up an experiment where the number of 
particles be a well controled quantity and the interactions strong. Nowadays it seems possible 



to construct quantum dot rings [0,0 which could be an ideal tool for this purpose. In 
this regard we would like to mention that recently, a fractional periodicity of (fro/ 4= has been 
observed in Auln rings prepared by e-beam lithography [F2U . 
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FIGURES 



FIG. 1. Energy versus magnetic flux for 4 particles in a 10 sites ring. For U = 1 it is 
observed that a 0o/2 periodicity begins to appear in the energy curve. For greater values of 
U, the system tends to a o /4 periodicity. 

FIG. 2. Energy versus magnetic flux for 5 particles in a 10 sites ring. For U = 1 there is 
no peaks in the curve. Peaks appears for U ~ 50 and in the large U limit the system tends 
to a 0o/5 periodicity. 

FIG. 3. Ground state energy as a function of flux for 4 particles in a 10 sites ring obtained 
with the U — oo Hubbard model. With open squares the energy in the F = 2 subspace, 
and filled squares for the energy in the F = 4 subspace. 

FIG. 4. Energy vs magnetic flux for 8 particles in a 12 sites ring obtained with the 
t-J z model for three values of J z . From up to down figures correspond to J z = 0,0.01,0.1. 
In each flux interval between sharp peaks, the ground state correspond to the numbers F 
indicated in the graph. 

FIG. 5. Energy vs magnetic flux for 8 particles in a 12 sites ring obtained with the t-J 
model for four values of J. 
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